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Direct Observation of the Extended Molecular Atmosphere of o Get by 
Differential Spectral Imaging with an Adaptive Optics System 
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Abstract 

We present new measurements of the diameter of o Get (Mira) as a function of wavelength in the 
2.2 /^m atmospheric window using the adaptive optics system and the infrared camera and spectrograph 
mounted on the Subaru Telescope. We found that the angular size of the star at the wavelengths of CO 
and H2O absorption lines were up to twice as large as the continuum photosphere. This size difference 
is attributable to the optically thick CO and H2O molecular layers surrounding the photosphere. This 
measurement is the first direct differential spectroscopic imaging of stellar extension that resolves individual 
molecular lines with high spectral-resolution observations. This observation technique is extremely sensitive 
to differences in spatial profiles at different wavelengths; we show that a difference in diameter much smaller 
than the point spread function can be measured. 

Key words: infrared: general, instrumentation: adaptive optics, stars: AGB and post-AGB, tech- 
niques: spectroscopic 



1. Introduction 

Mira-typc stars arc important in studying galactic evo- 
lution because the gases and dust ejected from them con- 
stitute a major contribution to galactic heavy element en- 
richment. However, detailed study of the physical struc- 
ture of the atmosphere of Mira-type stars was not feasible 
until recently, because the angular size of these stars is 
smaller than the resolution achieved using secing-limited 
optical and infrared telescopes. Although their deep ab- 
sorption features are clear evidence of the presence of 
molecular gases such as H2O and CO, their spatial dis- 
tribution in the outer atmosphere is not well understood. 

Recent observations using satellites and ground-based 
interferometers have greatly improved our knowledge 
of the spatial distribution of molecular gases. High 
spatial-resolution observations using infrared interferome- 
ters have been used to study stellar diameters (Ridgway et 
al. 1992; Dyck et al. 1996; Dyck, van Belle, & Thompson 
1998). Mennesson et al. (2002) compared the stellar size 
of cool M giants at K' and L' and found that stellar di- 
ameters are larger at longer wavelengths. They attributed 
this to emissions from an extended, partially transparent 
shell, although the source of the opacity was not identi- 
fied. Ireland et al. (2004) showed that measured stellar 
diameters are larger at molecular absorption bands in the 
near infrared region. Perrin et al. (2004) reported that the 



diameter is larger at the 2.39 /im band due to water vapor 
and CO molecules, and that molecular layers extend out 
to 2.2 R„. 

Based on spectroscopic observations using the Infrared 
Space Observatory (ISO), Yamamura, de Jong, & Cami 
(1999) modeled the outer atmosphere of Mira, a represen- 
tative late M-type star, which has a "hot" H2O shell at 
2.0 i?* with a column density of 3 x 10^^ cm~^ and a 
temperature of 2000 K, and a "cool" H2O shell at 2.3 
characterized by a column density of 3 x 10^° cm"'^ and a 
temperature of 1400 K. Here i?* denotes the radius of the 
stellar photosphere. Matsuura et al. (2002) shows that 
the very extended water vapor layer is a common feature 
in the Mira variables. 

Working from near-infrared (NIR) to mid-infrared 
(MIR) intcrferometric observation data, Ohnaka (2004) 
presented a two-layer model consisting of hot and cool 
H2O layers for three Mira variables and concluded that 
they have a hot 1800-2000 K layer between 1.5 and 1.7 i?* 
with a column density of 1 to 5 x 10^^ cm~^, and a cool 
1200-1400 K layer between 2.2 and 2.5 i?* with a column 
density of 1 to 7 x 10^^ cm^^. Weiner (2004) explained 
the enhancement of stellar size from NIR to MIR using a 
2000 K water vapor shell model. 

Dynamic modeling of the pulsation of Mira variables 
shows that the atmosphere extends well outside the photo- 
sphere where temperature is low enough to form molecules 
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(Hofner et al. 1998). Attempts to identify the molecular 
species contributing to the extension of the photosphere 
using multiwavelength intcrferometry in the J, iJ, and K 
bands were conducted with the Infrared Optical Telescope 
Array (IOTA) (Weigelt et al. 2003). Observations with 
narrow band filters in the K band were carried out with 
the Palomar Testbed Interferometer (Thompson, Creech- 
Eakman, & van Belle 2002). However, the spectral resolu- 
tion was not sufficient to resolve molecular lines. Recently, 
the Keck interferometer (spectral resolution R = 230) was 
able to observe the CO and water vapor bands at 2.0- 
2.4 /xm(Eisner et al. 2007). While those lines were not 
clearly resolved at that spectral resolution, their obser- 
vations start to show the molecular features contributing 
to the extension of the photosphere. The exposure time 
for intcrferometry has to be very short, and the telescope 
aperture sizes are typically small. Therefore, it is diffi- 
cult to increase the spectral resolution of an interferom- 
eter while maintaining a sufficiently high signal-to-noise 
ratio. 

In this paper, we use an instrument featuring adaptive 
optics (AO) correction that is suited to studying extended 
structures around a compact source with a spectral resolu- 
tion sufficiently high that it can resolve individual molec- 
ular lines. The diffraction- limited point spread function 
(PSF) of 8-m class telescopes is 32-56 mas at 1.25-2.2 /im, 
while the largest stellar disk angular size spanned by the 
nearby red super-giant a Ori is 44.2 mas at 2.2 /im (Dyck 
et al. 1992). Therefore, the spatial resolution is insuf- 
ficient to measure its stellar diameter by direct imaging 
under variable seeing conditions. Nevertheless, we can 
derive the wavelength dependence of the stellar diame- 
ter by conducting spectroscopic observations with AO. 
Diffraction-limited spectroscopy provides one-dimensional 
imaging of the stellar size at hundreds of spectral resolu- 
tion elements recorded simultaneously. By comparing the 
spatial profiles at some line features to those at the imme- 
diately adjacent continuum wavelength, one can achieve 
very sensitive measurements of the stellar diameter be- 
cause the PSF variation with time is essentially common 
and eliminated. We used this differential spectral imaging 
technique to investigate the stellar size in and out of the 
molecular lines to observe the extended molecular atmo- 
sphere of the late-type star o Get directly. 

The observation procedure used in differential spectral 
imaging is essentially the same as that of classical disper- 
sion spectroscopy. Some issues, however, must be care- 
fully considered. First, the correction by AO must be 
sufficiently good to deliver a PSF with a clear diffraction- 
limited core. Second, observations should be carried out 
on low-humidity nights because the water vapor in the 
stellar atmosphere is one of our target molecules and the 
PSF depends on the refraction indices of the telluric atmo- 
sphere, which vary considerably at these strong absorption 
lines. In our experience, as long as the humidity is below 
10%, which is not uncommon at Mauna Kea, the effect 
of the wavelength dependency of the refraction index is 
negligible. 



2. Observation 

o Get (Mira, M5e-M9e, 332-day variability period; 
Kholopov 1998) is one of the brightest late-type M giants. 
The outer structure of the star has been studied using 
space-based spectroscopy (Yamamura, de Jong, & Gami 
1999) and intcrferometry at optical and IR wavelengths 
(Ridgway et al. 1992; Haniff et al. 1995; Woodruff et al. 
2004). We conducted spectroscopic imaging of the object 
using the 36-element curvature-based adaptive optics sys- 
tem of the Subaru Telescope to achieve diffraction-limited 
resolution at 2 micron (Takami et al. 2004; Gacsslcr et al. 
2002). 

The observations were made on UT 2001 December 22 
during the commissioning run of the AO system using 
the Infrared Gamera and Spectrograph (IRGS) (Tokunaga 
et al. 1998; Kobayashi et al. 2000). The night was ideal 
for the planned observation as the humidity was low and 
the seeing conditions were favorable. We chose the K 
band covering the absorption bands from two abundant 
molecules, GO at 2.3 /im and II2O at 2.45 /im. 

A medium-resolution grism in the IRGS with 
23 mas/pixel was used with a 0''10 slit to provide spec- 
tra at a resolving power of 1000. The position angle of 
the slit was set along the east-west direction. The ob- 
ject itself served as a wavefront sensing reference for the 
AO system. The spectrograms were recorded by chopping 
the AO tip-tilt mirror by 3" along the slit. To avoid de- 
tector saturation from the very bright central object, we 
used a neutral density filter, and set the integration time 
to 0.26 s for a single exposure. One hundred exposures 
were co-added before readout to storage. The total time 
from the start of the first exposure to the end of the last 
exposure was 26 s including read-out after each integra- 
tion, which is sufficiently long to average out the speckle 
noise. A comparison star of similar spectral type (SW Get, 
M7 III; Kholopov 1998) was also observed. SW Get is 60 
times fainter than o Get at K and is used as a reference 
point source, since it is located eight times farther away 
from us than o Get. An early type spectroscopic stan- 
dard star was observed through a similar air mass as the 
program star to remove the telluric absorption lines. The 
spectroscopic flat fields were obtained at the end of the 
night with a halogen lamp. 

3. Data Analysis and Discussion 

We subtracted pairs of frames that had the spectral im- 
ages positioned 3" apart from each other on the IR array 
detector to remove the sky emission and dark-count pat- 
terns. The spatial profile at each spectral element along 
the slit was fitted with a three-component Gaussian func- 
tion to represent the diffraction core, the shoulder shape 
of the first Airy ring, and the seeing halo. This is the 
typical PSF profile with partially corrected AO observa- 
tion; therefore, we have introduced those components to 
our fitting calculation. There are no fixed parameters in 
the fitting. The free parameters arc the intensity, spec- 
tral width, and the wavelength for all three components. 
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The full-width at half-maximum (FWHM) of the spatial 
extent was measured against the synthetic profile at each 
wavelength. The spectra are slightly tilted along the pixel 
row by as much as 1%. However, rectification using a 
sub-pixel shift was not performed to avoid resampling ef- 
fects on the spatial profile. The effect of non-rectification 
is negligible. If it has any effect, the measured FWHM 
should have a periodic pattern at about a 100-pixel inter- 
val. However, we did not observe such a pattern in the 
derived FWHM. 

The spectra were extracted using the IRAF^ aperture 
extraction package. The wavelength was calibrated by 
maximizing the cross-correlation of the spectra with the 
modeled telluric atmospheric absorption curve calculated 
using the program ATRAN (Lord 1992). 

Figure 1 shows the spectra and spatial profiles of o Get, 
the target star, and SW Get, the reference unresolved 
star. The spatial extent of o Get is slightly greater at 
the molecular absorption lines than at the adjacent con- 
tinuum wavelength. In contrast, no variation is seen in 
the spatial profile of SW Get, measured in the same way 
as that of o Get. In Figure 2, we see that the increase in 
the measured FWHM of o Get is well correlated with the 
H2O opacity feature in the 2.43-2.48 /im region and the 
GO bandhead in the 2.29-2.41 /im region. This is clear 
evidence for extension of the water vapor and GO atmo- 
spheres outside the photosphere of o Get. Note that the 
FWHM of the reference star SW Get, an unresolved point 
source, is uniform and is not correlated with the absorp- 
tion line features. It shows that a difference in diameter 
much smaller than the point spread function can be mea- 
sured with this differential spectral imaging technique. 

The stellar radius is smaller than the continuum at the 
blue end of the GO bands at 2.29 and 2.32 /im. These do 
not seem to be real features, but we do not have a clear 
explanation for their origin. One possible explanation is 
that they could be artifacts caused by the sharp drop in 
flux at these wavelengths. Note, however, that this ef- 
fect is very small in the SW Get data with similar GO 
absorption features (Figure 2(b)). 

The size of the extended gas is estimated from the cur- 
rent data as follows. The measured apparent stellar size is 
the convolution of the real size with the diffraction-limited 
PSF at 2.2 /im for an 8.2-m telescope. Deconvolving the 
profile to evaluate the real stellar size directly from AO 
observations is not straightforward, because the PSF of 
the AO observations is time varying. Therefore the de- 
convolved stellar profile from a comparison star is not 
usually correct. To overcome the time variation prob- 
lem, we compare the stellar size at the wavelength of the 
molecular lines to that at the continuum wavelength in 
the same observation of Mira. That is, we use the stellar 
profile at the continuum wavelength as the simultaneous 
PSF reference. By adopting the photosphere size derived 
from other interferometric observations in the continuum, 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, under a cooperative 
agreement with the National Science Foundation. 



we can derive the stellar size measured at the molecular 
lines. 

Although the center-to-limb variation (GLV) of the in- 
tensity across the stellar disk has a complex profile, de- 
pending on the stellar model used (Scholz 2001; Jacob & 
Scholz 2002), we assumed a Gaussian function for sim- 
plicity. The Gaussian GLV is sufficient to estimate the 
FWHM of the molecular layer because we do not have 
any higher spatial resolution information. The apparent 
spatial width is then a square sum of the PSF size and 
the angular size of the star. 

We need to know the size of the PSF as a function of 
the wavelength, 6'pgF(A), to calculate the stellar diame- 
ter. However, the PSF measured with the comparison 
star in this observation is not useful because it can differ 
from that obtained during observation of the object due to 
changes in the seeing conditions. Instead, we estimate the 
system PSF using the following procedure. We assume 
a stellar photosphere without any molecular layer to be 
the size given by 0,o- The photosphere defines the stellar 
diameter at all wavelengths where there are no molecu- 
lar lines. The derived photospheric stellar size, 0obs_fit(A), 
should have the following relationship with the system 
PSF. 

^obs_fit(^) = ^*o(-^) + ^PSf('^)- 

From Figure 2, we obtained 0obs_fit(A) by fitting the 
lower envelope of the observed FWHM with a polynomial 
function. The lower envelope is our approximation to the 
continuum level of the star spectrum. The PSF size is 
then calculated with 6'pgp ( A) = 6'obs_fit W~^loW- We as- 
sume that the intrinsic diameter of the stellar photosphere 
0,o(A) is constant in the relevant wavelength region. The 
diameter 6'h.o of the photosphere at 2.22 /im, 20.8 +/- 
0.28 mas, was provided by fitting the visibility data of 
interferometric observation by Perrin et al. (2004). Their 
observation (phase (j)=0.20) was made one month prior to 
ours ((/>=0.30). Since Perrin et al. (2004) found only a 
3% increase in the photospheric diameter from 0=0.01 to 
0=0.2, we assume the same photospheric diameter. Since 
we normalize to the photospheric diameter, the absolute 
value is not critical. Then, the stellar diameters including 
molecular lines are 

^*(A) = v'eL(A)-[eo\s.fit(A)-e^o]- 

Figure 3 presents the stellar size so derived as a function 
of wavelength, and shows that for o Get, the molecular 
layers of water vapor and GO extend to twice the photo- 
spheric radii. We have examined the error in deriving the 
enhancement of the molecular layer radius from the photo- 
sphere. The error sources are the combination of random 
and systematic error. The random error is estimated to be 
+/- 1%, which comes from the noise level of the radius in 
the wavelengths of no molecular features. For the system- 
atic error, we have identified that the moderate spectral 
resolution is a major source. The size of 0obs_fit(A) might 
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be an overestimate because we do not resolve the indi- 
vidual molecular lines by our spectral resolution of R = 
1000. The bottom level of the molecular opacity curve 
smoothed with this spectral resolution is not zero. This 
means that the apparent continuum size we have derived 
by lower envelope fitting gives a larger value than the ac- 
tual size. This would lead to an underestimate of the 
extent of the molecular atmosphere. We have estimated 
this effect for the lines between 2.45 and 2.48 fim. When 
we look at the model opacity curve at this wavelength, the 
real apparent continuum size will be about 80 mas, while 
we set the lower envelope at this wavelength to 82.5 mas 
in Figure 2. This changes the estimated stellar size from 
2.12 to 2.33 R*. The uncertainty of the radius due to this 
effect is about 10%. 

The lower spectral resolution observation also provides 
another underestimate of the extent of the molecular layer, 
because the peaks of the opacity are also smoothed out 
and do not reach the peak positions of the individual lines. 

The error of continuum diameter of Mira taken from 
Pcrrin et al. (2004) is a source of the error in calculating 
the molecular layer diameter. From the above formula, 
the error of molecular diameter is estimated to be the 
same order of error of the continuum diameter, which is, 
approximately 1%. 

We compared our result to previous studies. The H2O 
atmosphere of o Get measured using the ISO spectroscopic 
observations (Yamamura, de Jong, & Cami 1999) extends 
to 2.3 (0=0.99). Weiner (2004) modeled the gas layer 
from the spectroscopic and interferometric data and de- 
rived a diameter of 2.4 (0=0.81). Perrin et al. (2004) 
got 2.0 i?, at 2.39 fim from their optical interferometer 
observations with IOTA. These observations gave num- 
bers similar to our measurements. Ireland et al. (2004) 
reported a diameter of 1.5 i?* at 2.46 ^m using Keck aper- 
ture masking observations, a value that is significantly 
smaller than our result. However, their spectral resolu- 
tion was much lower than ours. 

Observations at higher spectral resolution should re- 
duce these uncertainties and provide improved measure- 
ment for this type of study. The signal-to-noise ratio of 
the Subaru Telescope is sufficiently high to make obser- 
vations at a resolution ten times greater with the IRCS 
using this method. 

4. Summary 

We observed o Get with adaptive optics in the if-band 
with diffraction-limited spatial resolution. Using the dif- 
ferential spectral imaging technique, we confirmed that 
the stellar size measured at the molecular absorption lines 
is significantly larger than that measured at the adjacent 
continuum wavelengths. 

The effective diameters of the II2O and GO atmospheres 
of o Get are about twice that of the photosphere ob- 
served in the continuum. This is the first direct spec- 
troscopic imaging of the extended molecular atmosphere 
of a late-type star with sufficient spectral resolution to re- 
solve those lines. The size of the atmosphere measured is 



consistent with infrared interferometric observations and 
the spectroscopic study of the star by ISO. 

We thank all the staff and crew of the Subaru Telescope 
and NAOJ for their valuable assistance in obtaining these 
data and continuous support for the Subaru AO and IRGS 
construction. Miwa Goto is supported by a Japan Society 
for the Promotion of Science fellowship. 
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Fig. 1. (al) Blow-up of the 
2.4 ^m. The arrows point to 
an adjacent continuum region, 
sured spatial profiles, rectangles and crosses, extracted at the 
absorption line and at the continuum region marked in (al), 
respectively. The spectra appears wider at the bright con- 
tinuum wavelength. However, if we normalize to the peak 
intensity, the stellar size extends more at the molecular line 
wavelength. The solid red and green curves are three-compo- 
nent Gaussian functions fitted to the respective data points. 
The spatial profile at the absorption line is wider than that 
at the continuum region, (bl), (b2) Same as (al) and (a2), 
but for the PSF reference star SW Cet, which is eight times 
farther than o Cet and can be regarded as an unresolved point 
source. 
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Fig. 2. (a) The FWHM of the spatial profile as a function of 
wavelength for o Cet is presented in the upper panel and the 
intensity spectrum is shown in the lower panel. The FWHM is 
larger by up to 10% at the molecular lines relative to the con- 
tinuum wavelengths. The enhancement in the FWHM closely 
follows the molecular absorption features of H2O and CO. For 
comparison, the pure H2O opacity at 2000 K computed with 
the HITRAN/HITEMP database (arbitrary units) is given 
at the bottom of the upper panel (Rothmann et al. 1998). 
There is good correlation between the opacity plot and the 
measured stellar size at strong water vapor absorption bands 
longer than 2.42 fira. The spectrum of the opacity plot is 
smoothed with the spectral resolution R = 1000 of our ob- 
servation. We have fitted the lower envelope of the observed 
FWHM using a polynomial function ^?obs_flt(-^)- We regard 
the lower envelope as the apparent stellar size of the pho- 
tosphere without molecular lines convolved with the system 
PSF. (b) Same as (a), but for the comparison star SW Cet. 
Note that no FWHM enhancement is visible for this star. 
The uncertainty of the measurement is estimated from the 
variation in (b) to be approximately +/- 1 mas. 
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Fig. 3. Observed radius relative to tlie piiotospheric radius 
of o Get. The radius of the star extends up to twice the size 
of the photosphere at the H2O and CO hnes. The estimation 
is a lower hmit of the extension. The error of the radius 
from the overlapping molecular lines is estimated to be 10%. 
The effect of the low spectral resolution observation will be 
confirmed by observation with higher spectral resolution. 



